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ABSTRACT 

We present the results of a joint spectral analysis of RXTE PCA, ASCA SIS, and 
ROSAT PSPC data of the supernova remnant SN 1006. This work represents the first 
attempt to model both the thermal and nonthermal X-ray emission over the entire X-ray 
energy band. The thermal flux is described by a nonequilibrium ionization model with 
an electron temperature kT e = 0.6 keV, an ionization timescale n$t = 9 x 10 9 cm -3 s, 
^ ■ and a relative elemental abundance of silicon that is 10-18 times larger than the solar 

abundance. The nonthermal X-ray spectrum is described by a broken power law model 
with low- and high-energy photon indices T\ = 2.1 and 1^ = 3.0, respectively. Since the 
nonthermal X-ray spectrum steepens with increasing energy, the results of the present 
analysis corroborate previous claims that the nonthermal X-ray emission is produced 
by synchrotron radiation. We argue that the magnetic field strength is significantly 
larger than previous estimates of about 10 /iG and arbitrarily use a value of 40 //G to 
estimate the parameters of the cosmic-ray electron, proton, and helium spectra of the 
remnant. The results for the ratio of the number densities of protons and electrons 
(R = 160 at 1 GeV), the total energy in cosmic rays (E cr = 1 x 10 50 ergs), and the 
spectral index of the electrons at 1 GeV (r e = 2.14 ± 0.12) are consistent with the 
hypothesis that Galactic cosmic rays are accelerated predominantly in the shocks of 
supernova remnants. Yet, the remnant may or may not accelerate nuclei to energies as 
high as the energy of the "knee," depending on the reason why the maximum energy of 
the electrons is only 10 TeV. 



Subject headings: acceleration of particles — cosmic rays — ISM: individual (SN 1006) 
— radiation mechanisms: nonthermal — supernova remnants — X-rays: general 



1. Introduction 

The search for evidence of the origin of Galactic cosmic rays has been an active area of research 
for many decades. While little evidence exists about the sites at which very high energy nuclei are 
accelerated, the results of recent X-ray and gamma-ray observations indicate that at least some 
of the cosmic-ray electrons are accelerated in the shocks of supernova remnants (Koyama et al. 
1995, 1997; Allen et al. 1997; Tanimori et al. 1998; Allen, Gotthelf, & Petre 1999; Vink et al. 1999; 
Slane et al. 1999; Muraishi et al. 2000; Borkowski et al. 2000; Dyer et al. 2001a). For example, 
SN 1006 is one remnant for which there is evidence that cosmic-ray electrons have been accelerated 
to energies as high as about 100 TeV (Koyama et al. 1995; Tanimori et al. 1998). In this paper, 
measurements of the X-ray and radio emission of the remnant are used to determine the parameters 
of the nonthermal electron spectrum. Although there is no evidence to indicate that cosmic-ray 
nuclei are accelerated in SN 1006, relativistic electrons and nuclei are expected to be accelerated 
in a similar manner (Ellison & Reynolds 1991). Therefore, we estimate the parameters of the 
proton and helium spectra of the remnant. The results of this analysis show that the remnant is a 
significant source of Galactic cosmic rays (at least cosmic-ray electrons) and provide some support 
for the idea that Galactic cosmic rays are accelerated predominantly in the shocks of supernova 
remnants. 

Section 2 contains a description of some of the results of previous analyses of the X-ray data 
of SN 1006. The data and analysis techniques used for the present work are described in section 3. 
Section 4 includes a discussion of the results of the present analysis. The conclusions are reviewed 
in section 5. 



2. Observational X-Ray Record 

The first detection of X-ray emission from SN 1006 was reported by Winkler & Laird (1976). 
The results of analyses of X-ray spectral data of the remnant have been reported by many authors 
(Winkler & Laird 1976; Zarnecki & Bibbo 1979; Winkler et al. 1979; Toor 1980; Becker et al. 
1980; Pye et al. 1981; Galas, Venkatesan, &: Garmire 1982; Vartanian, Lum, & Ku 1985; Hamilton, 
Sarazin, & Szymkowiak 1986; Koyama et al. 1987, 1995; Leahy, Nousek, & Hamilton 1991; Ozaki 
et al. 1994; Willingale et al. 1996; Reynolds 1996; Laming 1998; Vink et al. 2000; Dyer et al. 2001a). 
As early as 1979, it was realized that the global X-ray spectrum of SN 1006 cannot be fitted with 
a single thermal or nonthermal model (Winkler et al. 1979). This realization is not surprising in 
retrospect because observations of the remnant with imaging X-ray devices show that the X-ray 
spectrum varies as a function of position in the remnant. The high-energy X-ray image of the 



remnant, like the radio image, is dominated by emission from the northeastern and southwestern 
rims (Pye et al. 1981; Willingale et al. 1996; Winkler & Long 1997). Since the spectra of the rims 
are more or less featureless, the emission from these regions is reported to be nonthermal (Koyama 
et al. 1995; Reynolds 1996; Dyer et al. 2001a). At energies of less than 0.5 keV, a significant amount 
of emission is observed from the interior of the remnant, especially in the southwest (Pye et al. 
1981; Willingale et al. 1996). Since the X-ray spectrum of the interior exhibits atomic emission-line 
features, a substantial amount of the emission from this region must be thermal bremsstrahlung. 

Several of the spectral analyses of SN 1006 included two spectral components. Invariably, one 
of the spectral components is thermal emission from an optically thin plasma. Such a model is 
needed to describe the emission-line features of oxygen, neon, magnesium, silicon, and perhaps 
iron (Koyama et al. 1995; Vink et al. 2000; Dyer et al. 2001a). Some analyses have included an 
additional thermal component, while others have included a nonthermal component. For example, 
Hamilton et al. (1986) modeled the spectrum using two thermal components. One component 
is used to describe the forward-shocked interstellar medium, and the other describes the reverse- 
shocked ejecta. To explain the nearly featureless nature of the global X-ray spectrum of the remnant, 
they used a model for the shocked ejecta in which most of the emission is produced by a nearly 
pure layer of fully ionized carbon. The elements more massive than carbon are assumed to be 
significantly underionized in the ejecta, to suppress the X-ray line emission that would otherwise 
be produced by these elements. However, these assumptions are inconsistent with the results of 
subsequent analyses. The results of an analysis of the ultraviolet absorption- line features along 
the line of sight to the S-M star suggest that most of the ejected silicon has already been shocked 
(Hamilton et al. 1997), and therefore may not be in a low ionization state. This result is supported 
by X-ray observations that reveal line emission from several highly ionized elements: oxygen, neon, 
magnesium, silicon, and perhaps iron (Koyama et al. 1995; Vink et al. 2000; Dyer et al. 2001a). 
Since the relative abundance of the shocked silicon producing the X-ray emission is much larger 
than the relative abundance of silicon in the solar system, most of the shocked silicon is ejecta. 
Furthermore, the results of the analysis of Laming (1998) suggest that the total amount of reverse- 
shocked carbon required by the model of Hamilton et al. may be unreasonably large for the remnant 
of a type la supernova (i.e., may be well in excess of the Chandrasekhar mass). For these reasons, 
the thermal model of Hamilton et al. provides a poor description of the X-ray emission of SN 1006. 

Willingale et al. (1996) present a nonthermal interpretation of the high-energy X-ray emission. 
They argue that the X-ray emission from the rims is produced by relativistic electrons that are 
beamed from an unobserved central compact object. This conclusion seems to be based largely on 
the X-ray morphology of SN 1006, because no evidence of a compact object or of particle jets in the 
remnant has been reported. For comparison, the pulsar and the pulsar jet of the Vela supernova 
remnant are observed to produce radio, X-ray, and gamma-ray emission (Markwardt & Ogelman 
1995; Yoshikoshi et al. 1997; Bock, Turtle, & Green 1998). Willingale et al. (1996) set an upper 
limit on the flux of a point source of 6.3 x 10~ 14 erg cm' 2 s -1 . The corresponding constraint on 
the blackbody temperature of a 10 km neutron star is Tbb < 6.9 x 10 5 K (at a distance of 1.7 kpc). 



Some of the models of cooling neutron stars indicate that the temperature of a neutron star is 
higher than this upper limit at an age of 990 yr (Slane &; Lloyd 1995). Furthermore, the sketchy 
record of the optical light curve of the supernova (Schaefer 1996), the lack of evidence of an OB 
association near SN 1006 (Hamilton et al. 1986), the large distance of the remnant from the plane of 
the Galaxy [z = +430(d/1.7 kpc) pc], measurements of the average expansion rate of the remnant 
(Moffett, Goss, & Reynolds 1993), and measurements of the relative abundances of the elements 
producing the ultraviolet absorption lines (Fesen et al. 1988) suggest that SN 1006 was produced by 
a type la supernova. Type la supernovae are generally not expected to produce compact objects. 
For these reasons, it seems unlikely that SN 1006 contains an undetected compact object that is 
responsible for the high-energy X-ray emission from the rims. 

Koyama et al. (1995) advance an alternative nonthermal explanation of the high-energy X-ray 
emission. They conclude that the emission from the rims is synchrotron radiation from electrons 
that have been accelerated to energies of about 100 TeV. Reynolds (1996) argues convincingly that 
this conclusion is the only plausible one and shows that the nonthermal X-ray emission of the 
rims can be described by a model in which both the radio and X-ray emission are produced by a 
common population of shock-accelerated electrons. The strong correlation of the features of a high- 
energy X-ray image and a radio image of the northeastern rim of SN 1006 (Winkler & Long 1997) 
supports this result. Furthermore, the detection of inverse-Compton TeV gamma-ray emission 
from the northeastern rim of the remnant confirms the presence of very high energy electrons 
in SN 1006 (Tanimori et al. 1998). Since the TeV photons are produced by electrons that have 
energies similar to the energies needed to produce X-ray synchrotron radiation and since the radio, 
X-ray, and gamma-ray spectral data can be fitted with simple models of the spectrum of shock- 
accelerated electrons (Mastichiadis &: de Jager 1996; Aharonian & Atoyan 1999), the conclusion 
that the nonthermal X-ray emission is produced by synchrotron radiation seems unavoidable. 



3. Data and Analysis 

To accurately model the shape of the entire X-ray spectrum of SN 1006, from 0.12 to 17 keV, 
we analyzed data obtained using the Proportional Counter Array (PCA) of the Rossi X-Ray Timing 
Explorer (RXTE) satellite, the Solid-state Imaging Spectrometers (SISs) of the Advanced Satellite 
for Cosmology and Astrophysics (ASCA), and the Position Sensitive Proportional Counter (PSPCB, 
hereafter just PSPC) of the Rontgen satellite (ROSAT). Between 1996 February 18 and 1996 
February 20, SN 1006 was observed for 21 ks at a location on the northeastern rim (02000 = 
15 h 4™0, c>2ooo = — 41°48') and for 20 ks at a location on the southwestern rim (02000 = 15 h l™8, 
^2000 = — 42°6') using the PCA. The PCA is a spectrophotometer that is comprised of an array of 
five coaligned proportional counter units that are mechanically collimated to have a field of view 
of 1° FWHM (Jahoda et al. 1996). The array is sensitive to photons that have energies between 
about 2 and 60 keV and has an energy resolution AE/E = 0.18 at 6 keV. The maximum on-axis 
collecting area is about 7000 cm 2 . The PCA data were screened to remove the time intervals 



during which (1) one or more of the five proportional counter units is off, (2) SN 1006 is less 
than 10° above the limb of the Earth, (3) the background model is not well defined, and (4) the 
pointing direction of the detectors is more than 0?02 from the nominal pointing direction in either 
right ascension or declination. The set of X-ray events that satisfies these four selection criteria 
includes 7 ks of the data of the northeastern pointing and 11 ks of the data of the southwestern 
pointing. Spectra were constructed for each of the two PCA pointings. The two spectra are quite 
similar because the field of view of the PCA is large enough to include the entire remnant at 
both pointing locations. Therefore, the PCA data are insensitive to modest changes in the shape 
of the spectrum of the remnant as a function of position in the remnant. The PCA background 
spectrum for SN 1006 was estimated using version 1.5 of the FTOOL pcabackest 1 . This version of 
pcabackest includes an estimate of the charged-particle and diffuse cosmic X-ray backgrounds based 
on the "VLE" count rate during observations of "source-free" regions of the sky and an estimate 
of the background associated with the decay of radioactive material that is activated when the 
spacecraft passes through the South Atlantic Anomaly. The tool has been found to produce slightly 
inaccurate estimates of the background associated with astrophysical sources. For example, when 
the background spectrum is subtracted from the spectral data of SN 1006, the resulting "source" 
spectrum of the remnant has a negative count rate at high energies. To compensate for this problem, 
the spectrum of the background associated with the South Atlantic Anomaly is multiplied by a 
factor of 0.2. This factor is determined by requiring the average difference between the source 
spectrum and the source model to be zero at energies between 27 and 71 keV. In this energy band, 
the total number of events associated with SN 1006 is expected to be negligibly small compared 
to the number of events associated with the background. The adjustment reduces the count rate 
of the 2-10 keV activation background from 1.0 to 0.2 counts s _1 for the top layer of anodes in 
the PCA. After applying this adjustment, the average 2-10 keV background count-rate for the top 
layer is estimated to be 17.2 ± 0.1 counts s _1 for both pointings. 

On 1993 August 19-20, 1993 September 13-15, and 1996 February 20-21, the instruments on 
ASCA (Tanaka, Inoue, & Holt 1994) were used to observe SN 1006 for 43 ks at a location in 
the center of the remnant (02000 = 15 h 2 m 48 s , ^2000 = — 41°55'46"), for 38 ks at a location on 
the northeastern rim (02000 = 15 h 3 m 32 s , ^2000 = — 41°46'25"), and for 57 ks at a location on the 
southwestern rim (02000 = 15 h 2 m 34 s , #2000 = — 42°2'58"), respectively. The SIS detectors are CCD 
devices (Burke et al. 1994) sensitive to photons that have energies between about 0.5 and 10 keV. 
The devices have an energy resolution AE/E = 0.06-0.08 (FWHM) at 1.5 keV and a field of view 
of 22' x 22'. The spatial resolution of the telescope < 1.5' (50% encircled-photon radius) for the SIS 
data used here. The maximum on-axis effective area of the system is about 220 cm 2 . The SIS data 
were screened using the standard "REV2" criteria that exclude time intervals (1) associated with 
passages through the South Atlantic Anomaly and regions of relatively high particles fluxes, (2) 
during which the bright limb of the Earth is in the field of view, and (3) during which the source 
is occulted by the Earth. The set of X-ray events that satisfies these selection criteria includes 23, 



For more information about FTOOLS, see http://heasarc.gsfc.nasa.gov/docs/software/ftools/ftools_menu.html. 



20, and 43 ks of data for the center, northeastern rim, and southwestern rim, respectively. The SIS 
spectra for these three spatially separate regions, which correspond to the central, northeastern, 
and southwestern ovals of Figure 1, have event rates of 0.21, 0.89, and 0.66 counts s _1 , respectively. 

On 1993 January 28-29, SN 1006 was observed for 6 ks at a location near the center of the 
remnant (a 20 oo = 15 h 2 m 55 s , 5 2 ooo = -41°55'12") using the PSPC. The PSPC is a multiple-wire 
proportional counter sensitive to photons that have energies between about 0.1 and 2 keV. The 

o 

energy resolution AE/E = 0.43 at 0.93 keV. The telescope has a field of view of 2 . The spatial 
resolution < 30" (50% encircled-photon radius) for the PSPC data used here. The maximum on-axis 
effective area of the system is about 260 cm 2 (Pfeffermann et al. 1987). Spectra were constructed 
using the PSPC data for three spatially separate regions: the bright northeastern rim, the bright 
southwestern rim, and the center of the remnant (the three rectangular regions of Fig. 1). The 
total PSPC event rates for these three spectra are 5.6, 6.1, and 11.9 counts s _1 , respectively. The 
PSPC background for SN 1006 is negligibly small. Plucinsky et al. (1993) report that the event 
rate associated with unrejected cosmic-ray particles is 4 x 10~ 6 counts s" 1 keV -1 arcmin -2 in the 
pulse-height range 18 < PHA < 249. In this case, the cosmic-ray event rate over the entire remnant 
is only about 3 x 10~ 3 counts s _1 keV -1 . Since this rate is considerably smaller than the PSPC 
count-rate for SN 1006 at all energies from 0.12 to 2 keV, the backround is neglected. 

A joint fit to the RXTE PCA, ASCA SIS, and ROSAT PSPC data was performed using 
version 10.0 of the spectral-fitting software package XSPEC 2 . Since Koyama et al. (1995) and others 
report evidence of both thermal and nonthermal emission, the spectral fits were performed using 
several combinations of thermal and nonthermal emission components. The thermal components 
include a bremsstrahlung model (Brem), a model of a thin thermal plasma in collisional ionization 
equilibrium (RS; Raymond & Smith 1977), and a model of a thin thermal plasma that is not in 
ionization equilibrium (NEI; Hamilton, Sarazin, & Chevalier 1983). The nonthermal components 
include power law (PL) and broken power law (BPL) models. For each fit, only the absolute 
normalization of each spectral component was allowed to vary from one region to another. While 
the spectra most likely vary as a function of position in the remnant, we find no evidence of a 
significant variation in the shape of the nonthermal spectra of the different regions. Furthermore, 
Koyama et al. (1995) and Vink et al. (2000) report no evidence of spatial variation in the emission 
of the interior of the remnant. 

Table 1 summarizes the results of some of the fits. Note that the values of x 2 l v > 1 f°r an 
of the fits. This result suggests that the value of x 2 may be sensitive to systematic errors in the 
responses of the detectors as well as differences between the spectral models and the spectral data. 
Therefore, the x 2 test statistic does not provide an unambiguous measure of the goodness of fit of 
the models. However, if the contributions of the systematic errors to % 2 do not change from one 
fit to another, the differences in x 2 provide a useful measure of the relative goodness of fit of the 
different models. In this case, Ax 2 = X 2 ~ Xmin ^ s distributed as a x 2 variable with the number of 
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degrees of freedom equal to the number of free parameters in the fit (Lampton, Margon, & Bowyer 
1976). For example, 26 free parameters are used for the best-fit model of Table 1. Therefore, the 
values of Ax 2 = 28.9, 35.6, 39.3, and 50.6 are expected to correspond to the 1 a, 90%, 2 a, and 
3 a confidence levels, respectively. These values of Ax 2 are used in this paper to determine the 
confidence intervals of the fit and inferred parameters. Since the differences between the value of 
X 2 for the BPL+NEI model and the values of x 2 f° r the other sets of models > 93, the BPL+NEI 
model seems to provide a significantly better fit to the entire broadband spectral data than the 
other sets of models listed in Table 1. Table 2 lists the fit parameters of the BPL+NEI model. 
The spectral data and the best-fit model are shown in Figures 2-4. In general, the model fits the 
0.12-17 keV data quite well. 

The thermal component of the best-fit model is a nonequilibrium ionization model of a thin 
thermal plasma that has different electron and ion temperatures (Hamilton et al. 1983). The shape 
of the X-ray continuum is specified by the temperature associated with the forward shock, Tf s , and 
the parameter n^Eo, where no is the ambient density of interstellar hydrogen and Eq is the initial 
kinetic energy of the ejecta. This model is implemented in XSPEC version 10.0 using a fixed grid 
in the T^-UqEq plane. The temperature is quantized by factors of 10 ' 25 from 10 6 ' 25 to 10 8,25 K 
and UqEq is quantized by factors of 10 from 10 48 to 10 53 ergs cm" 6 . The best-fit values of these 
two parameters are Tf s = 10 80 K and UqEq = 10 50 ergs cm" 6 (Fig. 5). Figure 5 includes curves 
along which the ionization timescale, not, and the characteristic electron temperature, kT e , are 
constant. The elongated confidence level contours in this figure nearly lie along curves of constant 
ionization timescale. The best-fit values of Tf s and UqEq are driven by the thermal emission of 
SN 1006 because the same values for these two parameters (within the statistical uncertainties) are 
obtained if the same set of spectral models is fitted to only the ASCA SIS spectrum of the center 
of the remnant. The relative elemental abundances of the elements oxygen, neon, magnesium, and 
silicon are included as free parameters in the fit because the ASCA SIS spectrum of the central 
region of SN 1006 exhibits atomic emission-line features associated with these elements (Koyama 
et al. 1995). The abundances of sulphur and calcium are fixed to be the same as the abundance of 
silicon because this situation is more or less consistent with the type la supernova model of Nomoto, 
Thielemann, Sz Yokoi (1984). Since the value of x 2 1S n °t significantly reduced if the abundances 
of helium, carbon, nitrogen, iron, or nickel are included as free parameters in the fit, these five 
abundances are set to be the same as the relative abundances of these elements in the solar system 
(Anders & Grevesse 1989). The fitted abundances of oxygen and magnesium are consistent with 
the solar abundances of these elements. The abundance of neon in SN 1006 is less than 80% of the 
relative solar abundance of neon. The fitted relative abundance of silicon (and sulphur and calcium) 
is 10-18 times larger than the solar abundance of silicon. The uncertainties of the abundances are 
somewhat larger than the statistical uncertainties specified in Table 2 because (Hamilton et al. 
1983) note that the intensities used in the NEI model for the strong emission lines are generally 
uncertain by at least a factor of 2. 

A broken power law model is used here to approximate a gradually steepening nonthermal X- 



ray spectrum. Such a spectrum is expected if the nonthermal emission is produced by synchrotron 
radiation (Reynolds 1996). As shown in Figure 6, the nonthermal spectrum of SN 1006 steepens 
significantly with increasing energy. The difference between the high-energy and low-energy photon 
indices of the broken power law > 0.7 at the 3 a confidence level. The value of the break energy of 
the broken power law (E(, = 1.85±0.18 keV) is not physically meaningful. It does not correspond to 
a feature in the nonthermal X-ray spectrum of SN 1006. This value is determined by the transition 
from the relatively high count rate portion of the ROSAT PSPC spectrum to the relatively high 
count rate portion of the RXTE PCA spectrum. 



4. Discussion 

4.1. Thermal Emission 

Only one thermal emission model is used in the present analysis to describe the properties of 
the X-ray-emitting material in SN 1006. The use of one thermal component is undoubtedly an 
oversimplification, because both the forward-shocked interstellar material and the reverse-shocked 
ejecta produce X rays. The large fitted abundance of silicon suggests that most of the shocked 
silicon is ejecta. Yet the inferred amount of swept-up mass implies that most of the shocked mass 
is interstellar. For these reasons, the use of a single thermal model to describe the global thermal 
emission properties of the remnant may provide a poor description of any given volume element 
of the X-ray-emitting material. Nevertheless, the use of a single thermal model is sufficient to 
accurately characterize the properties of the nonthermal X-ray emission (our primary goal) , because 
the X-ray emission of the remnant is dominated by nonthermal emission at energies above 1 keV, 
because the results for the broken power law component are insensitive to the details of the single 
thermal component used in the fits (compare the results obtained for the broken power laws of 
the "BPL+NEI" and "BPL+RS" models listed in Table 1), and because Koyama et al. (1995) and 
Vink et al. (2000) report no evidence of spatial variation of the thermal X-ray spectrum of the 
interior of SN 1006. 

Models that include two thermal components and one nonthermal component were fitted to 
the data. However, the results for the fit parameters of these three-component models showed that 
the parameters of one of the two thermal components are essentially unconstrained. 

Vink et al. (2000) favor a model for the X-ray emission of the interior of SN 1006 that includes 
two thermal components instead of one thermal component and one nonthermal component. They 
find that the temperatures of the hotter thermal components kT e > 3 keV for both the northern 
and southern regions of the interior. These temperatures are quite high for the observed Galactic 
supernova remnants. Of the young, shell-type remnants, only Cassiopeia A is reported to have a 
thermal plasma with a temperature this high. Hamilton & Sarazin (1984) show that it may be 
possible for a reverse-shocked silicon-dominated plasma to have electron temperatures this high, 
but such a metal-rich plasma would cool quickly. Therefore, the cooling time would have to be 



comparable to the age of the remnant or longer. Otherwise, the reverse-shocked electrons may 
be too cool to produce much X-ray emission. The model of Vink et al. (2000) cannot explain 
nonthermal X-ray emission from the interior of the remnant. Since Dyer et al. (2001b) show that 
about 25% of the radio emission at 843 MHz is produced in the interior of SN 1006, a similar fraction 
of nonthermal X-ray synchrotron emission is expected from this region. Our fits to the X-ray data 
with a model that includes one thermal and one nonthermal component suggest that 13 ± 7% of 
the nonthermal X-ray emission is produced in the region. Since this fraction is comparable to the 
fraction of the radio emission from the interior of SN 1006, we favor a model for the interior (and 
the rims) that includes one thermal component and one nonthermal component instead of two 
thermal components. It should be emphasized that the principal results of the present analysis 
concerning the properties of the cosmic rays in SN 1006 are insensitive to this choice, since only a 
small fraction of the high-energy X-ray emission is produced in the interior of the remnant. 

The parameters of the best-fit model are listed in Table 2. Table 3 includes a list of the inferred 
values of the characteristic energy of the thermal electrons, kT e (see Fig. 5), the velocity of the 
forward shock, Vf s , the radius of the forward shock, rf s , the distance of the remnant, d, the ionization 
timescale, not (see Fig. 5), the age of the remnant, t, the density of interstellar hydrogen, no, and 
the swept-up mass, M s (Hamilton et al. 1983). These inferences are based on the assumptions that 
SN 1006 is in a Sedov phase, that the initial kinetic energy of the ejecta is 10 51 ergs, that the X-ray 
emission is dominated by forward-shocked interstellar material, and that the distribution of the 
ambient material around SN 1006 is homogeneous. 

The best-fit value for the absorption column density of interstellar hydrogen nu = (5.6 ±0.6) x 
10 20 atoms cm~ 2 (at the 90% confidence level). This value is consistent with the result of Schaefer 
(1996) [njj = (6.1 ± 0.6) x 10 20 atoms cm' 2 ], who analyzed a compilation of results from several 
sources. The relative elemental abundances of oxygen, neon, magnesium, and silicon are consistent 
with the results of Koyama et al. (1995), but Vink et al. (2000) obtain a somewhat lower abundance 
for oxygen. The inferred energy of the thermal electrons (kT e = 0.58l ' 27 keV at the 90% confidence 
level) is similar to the results obtained by Vink et al. (2000; kT e = 0.71 ± 0.15 and 0.78 ± 0.09 keV) 
and Dyer et al. (2001a; kT e = 0.60l ' 52 ), but it is somewhat lower than the result of Koyama 
et al. (1995; kT e = 1.6 keV). The results of an analysis of the ultraviolet emission-line data of 
the northwestern rim favor a low electron temperature (kT e < 0.05kT p = 0.05-^.m p vf s = 0.7 keV), 
but the results cannot exclude temperatures as high as kT e = 2.9 keV (Laming et al. 1996). The 
inferred shock velocity (vf s = 2700^ 1700 km s _1 at the 90% confidence level) is consistent with the 
results of Kirshner, Winkler, & Chevalier (1987; V{ s = 2800-3870 km s' 1 ), Raymond, Blair, & Long 
(1995; V{ s = 2300 km s _1 ), and Laming et al. (1996; Vf s = 2600 ± 300 km s _1 ) for the northwestern 
rim of SN 1006, but Dwarkadas & Chevalier (1998) argue that the average forward shock velocity 
is 4000 ±500 kss- 1 . 

The inferred distance (d = 1.4^ 01 kpc at the 90% confidence level) is consistent with the result 
of Schaefer (1996; d = 1.59 ± 0.13 kpc) and with a distance estimate that is based on the velocity 
and proper motion of the forward shock in the northwest. As described above, Laming et al. 
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(1996) estimate that the velocity of the forward shock in the northwest Vf s = 2600 ± 300 km s -1 . 
Since neutral atoms have very short lifetimes after having passed through the shock into the hot 
postshock gas, the Ha filaments trace the location of the forward shock (Chevalier & Raymond 
1978; Long, Blair, & van den Bergh 1988; Winkler &: Long 1997). The results of analyses of 
the mean proper motion of the filaments along the northwestern rim of SN 1006 suggest that 
p, = / /30 ± / /04 yr -1 (Long et al. 1988) or p, = 0"39 ± 0"06 yr -1 (Hesser & van den Bergh 1981). 
The weighted average of these two results is 0'/33 ± / /03 yr -1 . Therefore the implied distance of 
SN 1006 d = Vf s /(p) = 1.66 ± 0.24 kpc. These distance estimates are consistent with the upper 
and lower limits on the distance (Hamilton et al. 1986, 1997; Wu et al. 1993; Burleigh et al. 2000). 
Collectively, the results suggest that 1.4 < d < 2.0 kpc. 

The ionization timescale (not = 240-650 cm -3 yr), age (t = 760-6300 yr), and ambient 
density of interstellar hydrogen (no = 0.10-0.37 cm -3 ), are consistent with the ionization timescale 
of Koyama et al. (1995; not = 300 cm -3 yr), the known age of the remnant, and the ambient 
density estimates of Toor (1980; no = 0.3 cm" 3 ) and Vink et al. (2000; no ~ 0.1 cm -3 ). Willingale 
et al. (1996) find a somewhat larger ambient density [no = 0.52 ± 0.08(d/1.7 kpc) 1 ' 2 cm -3 ], and 
Dwarkadas & Chevalier (1998) report a somewhat lower density (no = 0.05-0.1 cm -3 ). Although 
the remnant is interacting with a relatively dense environment in the northwest (no ~ 1 cm -3 ; 
Winkler & Long 1997), the average ambient denisty at a Galactic height of 430 pc is about 0.06 cm -3 
(Boulares & Cox 1990). Therefore, the average ambient density around SN 1006 may be about 
0.1 cm -3 . 

The inferred amount of shocked material M s = 8-67 Mq. However, this range of values is 
an overestimate because the relative abundance of helium used in the best-fit model (Table 2) is 
the same as the relative abundance of heilium in the solar system. Since Savedoff &: Van Horn 
(1982) and Burleigh et al. (2000) report much smaller densities of helium in the neighborhood of 
SN 1006, a more appropriate abundance may be about 0.2 times the solar abundance. In this 
case, the best-fit amount of shocked material M s = 6.8 Mq, and the 90% confidence level range is 
6-47 Mq. This range is consistent with the mass estimates of Pye et al. (1981; M s = 5-15 Mq) 
and Vink et al. (2000; M s = 8.3 ± 0.8 Mq). However, the results of Dwarkadas & Chevalier (1998) 
suggest a somewhat smaller mass (M s = 3-5) and Willingale et al. (1996) obtain a much smaller 
mass [M s = l(d/1.7 kpc) 1 / 2 Mq]. With the exeption of the estimate of Willingale et al. (1996), 
the mass estimates are substantially larger than the Chandrasekhar mass. Therefore, most of the 
hot shocked material is probably forward-shocked interstellar matter. 

Unlike the parameters inferred from the thermal component of the BPL+NEI model, most of 
the parameters inferred from the thermal component of the BPL+RS model are inconsistent with 
the results of other observations. For example, the inferred velocity of the forward shock Vf s = 
(16/cT/3 / u) 1 ' 2 = 395-420 km s -1 and the inferred ambient density no = 0.006-0.014 atoms cm -3 . 
These estimates are about an order of magnitude smaller than the results of previously published 
analyses. The values inferred for most of the other parameters listed in Table 3 are also unrealistic 
for the Raymond-Smith component of the BPL+RS model. Therefore, the Raymond-Smith model 
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provides a poor description of the thermal emission of SN 1006. The differences in the physical 
properties of the Raymond-Smith and nonequilibrium ionization models in XSPEC are that the 
Raymond-Smith model is based on the assumptions that the plasma is in a state of collisional 
ionization equilibrium and that the electron and ion temperatures have equilibrated. Neither of 
these two conditions is assumed for the nonequilibrium ionization model. In fact, at least one 
of the assumptions is not appropriate for SN 1006. Laming et al. (1996) show that the electron 
temperature is significantly smaller than the temperature of the ions in the northwest (kT e < 
0.2kTi). This difference in the temperatures could explain why the parameters inferred from the 
Raymond-Smith model are inconsistent with other results. 

Table 4 lists estimates of the kinetic energy of the matter in the remnant, E^ n , the thermal 
energy of the shock-heated protons, EkT p , and electrons, EkT e , the energy in the magnetic field, 
Eb, and the energy in cosmic rays, E cr . Some of these estimates are based on the assumptions 
that the remnant is in a Sedov phase (Hamilton et al. 1983) and that the initial kinetic energy of 
the ejecta Eq = 10 51 ergs. The scaling of the estimates with Eq is indicated in the table. However, 
the results for SN 1006 may be insensitive to the assumption that the remnant is in a Sedov phase. 
For example, Dwarkadas & Chevalier (1998) suggest that SN 1006 has not yet reached the Sedov 
phase. The results of their analysis suggest that M s pss 4 M Q , no ~ 0.07 cm" 3 , rf s = 8.5 pc, 
and t>f s = 4000 km s -1 . If these values (instead of the values in Table 3) are used to compute 
the parameters in Table 4, E\^ lTL = 4 x 10 50 ergs and E^x = 4 x 10 50 ergs. Therefore, the same 
distribution of energy is obtained using the results of Dwarkadas & Chevalier (1998) because the 
smaller estimates of the amount of shocked material and the ambient density are offset by the larger 
estimates of the shock radius and velocity. 



4.2. Nonthermal Emission 

The high-energy X-ray emission from the northeastern and southwestern rims of SN 1006 is al- 
most certainly nonthermal. An attempt to describe this emission using only thermal bremsstrahlung 
emission (Hamilton et al. 1983) was unsuccessful (Laming 1998). As described in section 2, the 
suggestion that the emission is associated with electrons beamed from an unseen central object is 
untenable. Reynolds (1996) argues convincingly that the only plausible explanation of the emission 
is synchrotron radiation from electrons accelerated to energies of about 10-100 TeV. This conclu- 
sion is supported by the results of the present analysis, which show that the nonthermal spectrum 
steepens with increasing energy (Fig. 6). 

Figure 7 displays a compilation of the radio flux-density results of Kundu (1970), Milne (1971), 
Milne & Dickel (1975), Stephenson, Clark, & Crawford (1977), and Roger et al. (1988). The radio 
data can be fitted with a power law spectrum: S v = (17. 9± l.l)(i//l GHz) -a57=l:0 ' 06 Jy. The results 
of analyses of data obtained using the IRAS (Arendt 1989), ROSAT PSPC, RXTE PCA, EGRET 
(Hartman et al. 1999), CANGAROO (Tanimori et al. 1998), and JANZOS (Allen et al. 1995) 
instruments are also shown. The figure includes estimates of the photon spectra associated with 
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synchrotron radiation, inverse Compton scattering of the cosmic microwave background radiation, 
nonthermal bremsstrahlung, and the decay of neutral pions Gaisser, Protheroe, & Stanev (1998). 
The synchrotron spectrum is extrapolated from the radio data to X-ray energies assuming that the 
electron spectrum is as described in section 4.3. The nonthermal bremsstrahlung and n° spectra 
were computed assuming the density of protons in the remnant is 0.7 cm -3 . The nonthermal X-ray 
spectrum is consistent with a model of synchrotron radiation, but it is not consistent with models 
of nonthermal bremsstrahlung or inverse Compton scattering. The TeV gamma-ray data can be 
described by inverse Compton scattering. 

As illustrated in Figure 7, the detection of X-ray synchrotron emission from SN 1006 requires 
the emission of synchrotron radiation at all energies between radio frequencies and X rays. Winkler 
& Long (1997) report the detection of faint diffuse emission (~ 1.0-2.5x 10 -17 erg cm' 2 s -1 arcsec -2 ) 
in the Ha band around the entire circumference of the remnant. Reynolds (1996) notes that the 
surface brightness of the synchrotron emission is about 28-32 mag arcsec -2 in the visible band. 
This range corresponds to a flux of 0.1-6 x 10 -17 erg cm -2 s -1 arcsec -2 in the Ha band. Therefore, 
some of the faint diffuse Ha flux may be produced by synchrotron radiation from electrons that 
have energies of about 1 TeV. The measurement of the optical (or infrared) synchrotron flux would 
be very important because it would corroborate the claim that the high-energy X-ray emission is 
synchrotron emission and could help determine the shape of the cosmic-ray electron spectrum of 
SN 1006. 

As shown in Figure 7, the synchrotron and inverse Compton spectra of SN 1006 have similar 
shapes. Gamma-ray emission is detected near the peak of the inverse Compton spectrum. The 
electrons that produce this emission are the same electrons that produce the emission near the 
peak of the synchrotron spectrum. Therefore, the gamma-ray and nonthermal X-ray emission are 
produced by electrons that have similar energies. For this reason, the nonthermal X-ray and gamma- 
ray images are expected to be similar. Tanimori et al. (1998) report the detection of gamma rays 
from only the northeastern rim of SN 1006. The 95% confidence level upper limit on the gamma-ray 
flux from the southwestern rim (1.1 x 10 -12 cm -2 s -1 at energies above 1.7 TeV) is no larger than 
about 40% of the flux of the northeastern rim [(4.6 ± 0.6 ± 1.4) x 10 -12 cm -2 s -1 at energies above 
1.7 TeV]. However, the nonthermal X-ray fluxes of the northeastern and southwestern rims are 
comparable. This discrepancy can be explained if the electron spectra, magnetic field strengths, 
or synchrotron and inverse Compton emission- angle distributions of the two rims are different. 
Of these possibilities, the X-ray data can be used to search for differences in the shapes of the 
electron spectra by searching for differences in the shapes of the synchrotron spectra of the two 
rims. Figure 8 shows a comparison of the nonthermal X-ray spectra of the two rims. For the 
purposes of this comparison, only the ROSAT PSPC data were used. The PSPC spectra of each 
rim were fitted with a model that includes a power law component (dN/dE = F(E/l keV) -r ) 
and the thermal component described in Table 2. The ratio of the nonthermal X-ray fluxes of the 
southwestern and northeastern rims f(E) = (Fsw/Fne)(E/1 keV) ne sw . The 1, 2, and 3 a 
confidence level contours for the parameter space defined by -Fsw/^ne and Tne — Tsw are shown 
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in Figure 8. Since the point denned by Fsw/^ne = 1 and Tne — Tgw = is well inside the 1 a 
confidence contour, there is no reason to believe that the shapes of the nonthermal X-ray spectra of 
the two rims are different in the energy range to which the ROSAT PSPC is sensitive. The results 
of the present analysis of the RXTE PC A data of the two rims also suggest that the shapes of the 
nonthermal X-ray spectra of the two rims are consistent with each other in the energy range to 
which the PC A is sensitive. Figure 8 includes curves along which the ratio / is a constant. Here, 
an energy E = 0.1 keV is used to compute the value of /, because we expect that the synchrotron 
emission at 0.1 keV and the TeV emission are produced by electrons that have the same energies. 
At the 3 a confidence level, the nonthermal flux from the southwestern rim at 0.1 keV is greater 
than 50% of the nonthermal flux from the northeastern rim at this energy. At face value, this 
lower limit seems to be inconsistent with the upper limit of 40% on the ratio of the gamma-ray 
fluxes. However, a significant difference between the electron spectra of the two rims produces a 
significant difference in the corresponding inverse Compton spectra, but only a modest difference 
in the synchrotron spectra. Therefore, the modest constraint on the ratio of the inverse Compton 
fluxes corresponds to a poor contraint on the on the ratio of the synchrotron fluxes of the two rims, 
and the results of the present analysis cannot exclude the possibility that the differences between 
the X-ray and gamma-ray images are due to differences between the electron spectra of the two 
rims. Perhaps this issue will be settled after the TeV and X-ray spectra of the two rims are better 
measured. 



4.3. Cosmic Rays 

The available radio and X-ray data for SN 1006 can be used to infer the parameters of the 
spectrum of the cosmic-ray electrons producing the synchrotron emission. For simplicity, the rela- 
tivistic electron spectrum is assumed to have the form dn e /dE = A e E e exp(—E/e e ). This form 
is the high-energy limit of Bell's (1978) expression: 

dn/dE = Ag(E,m,T,e) (1) 

= A(E + mc 2 )(E 2 + 2mc 2 £r (r + 1)/2 exp(-£/e), (2) 

where E [= (7 — l)mc 2 ] is the kinetic energy of the particles and the exponential cutoff is added to 
the formula of Bell. The electron spectral index T e can be determined from the spectral index of the 
radio data: T e = 2q + 1 = 2.14 ± 0.12. However, estimates of the normalization factor A e and the 
exponential cutoff energy e e of the electrons cannot be uniquely determined from the synchrotron 
data alone because the normalization (S u ex A e B 1+a ) and the roll-off energy (.E ro ii oc e 2 B) of the 
synchrotron spectrum depend on the strength of the magnetic field. (Note that E TO u = 0.1 keV 
in Fig. 7.) Since the normalization of the inverse Compton spectrum produced by the relativistic 
electrons depends on A e (and the known spectrum of the cosmic microwave background radiation), 
but not B, the TeV gamma-ray data can be used to determine A e . Then the flux and the roll- 
off energy of the synchrotron spectrum can be used to determine B and e e , respectively. For 
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example, if the cosmic-ray electrons and the magnetic field of SN 1006 have effective volume- 
filling factors f e = Jb = 0.25, if the volume of the remnant V = ^irO^d? = 5 x 10 58 cm 3 , and 
if the inverse Compton emission is dominated by scattering of the cosmic microwave background 
radiation, A e = 8.6 x 10~ 9 cm" 3 GeV r_1 , B = 10 //G, and e e = 20 TeV. Similar estimates of the 
magnetic field strength are reported by Tanimori et al. (1998; B = 6.5±2 (iG), Aharonian & Atoyan 
(1999; B = 10 /txG), and Dyer et al. (2001a, B = 10 fiG). However, these estimates of the magnetic 
field strength seem rather small. If the field strength at a Galactic height of 430 pc is similar to the 
strength in the plane (Boulares & Cox 1990), a field of 10 /iG can be little more than a compressed 
ambient field. Yet, the large and abrupt change in the radio surface brightness in the regions of the 
bright rims of SN 1006 (Reynolds & Gilmore 1986), the presence of a radially oriented component 
of the field (Reynolds & Gilmore 1993), and estimates of the total energy in cosmic-ray particles 
(described below) suggest that the field has been significantly amplified. Perhaps one or more of 
the assumptions used to determine the strength of the magnetic field is wrong. In particular, the 
effective volume-filling factor of the magnetic field may be smaller than the filling factor of the 
cosmic rays (Jun & Jones 1999). In this case, the magnetic field strength is underestimated. For 
example, if f e = 0.25 and f B = 0.1, A e = 2.4 x 10~ 9 cm~ 3 GeV r ~\ B = 40 fiG, and e e = 10 TeV. 
This value for the magnetic field lies comfortably between the compressed value (B ~ 10 fiG) and 
the value obtained using the minimum-energy condition {B ~ 100 //G). For the remainder of the 
discussion, it is assumed that B = 40 fiG and that the cosmic-ray electron spectrum is specified 
by equation (2), where A e , T e , and e e are listed in Table 5. The inferred spectrum of cosmic-ray 
electrons is shown in Figure 9. 

Although there is no evidence that cosmic-ray nuclei are accelerated in SN 1006 (cf. Aharo- 
nian & Atoyan 1999), Ellison & Reynolds (1991) suggest that relativistic electrons and nuclei are 
accelerated in a similar manner. Therefore, we estimate the parameters of the cosmic-ray nuclei 
spectra. For simplicity, only hydrogen and helium nuclei are considered. The nonthermal spectra of 
these particles are assumed to have the same functional form as equation (2). The spectral indices 
of the protons and alpha particles are assumed to be the same as the spectral index of the electrons 
(T p = Thc = r e = 2.14). The exponential cutoff energies of the different cosmic-ray particles are 
assumed to be related by the magnetic rigidity of the particles (i.e., e p = ^en e = e e = 10 TeV). 
This relationship is appropriate if the maximum energy of the electrons is limited by the escape 
of the particles from the remnant (Reynolds 1996; Dyer et al. 2001a), but it is not appropriate if 
e e is limited by radiative losses. The values of A p and ^He are determined using the value of A e 
and the technique described below. The ratio of the number density of protons to helium nuclei 
n p '■ n Hc = 1 : 0.02 in the neighborhood of SN 1006 (Savedoff & Van Horn 1982; Burleigh et al. 
2000). If the shocked protons and helium atoms are fully ionized, then n e = 1.04n p . (The electron 
contribution from other nuclei is expected to be small.) Therefore, n p : njjc ■ n e = 1 : 0.02 : 1.04. 
The fractional numbers of nonthermal protons and helium nuclei are assumed to be the same as 
the fractional number of nonthermal electrons (i.e., r\ = njj r /n p = n| r e /riHe = n c * /n e ). Using this 
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relationship and the equations 

rip = A p dEg(E,mp,T p ,e p ) (3) 

° J - y min,p 

and 
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nf = A e dEg(E,m e ,T e ,e e ) (4) 

-Cjmin.f 5 . 



yields 



J*™^dEg(E,m e ,r e ,e e ) 

A = A m "'' e (5) 

p e n e j^v dE{ET) 

° D mm,p c l l 

= 1.1 x 1(T 6 cm" 3 GeV r_1 , (6) 

where £ mi n,e = 0.58 keV, ^ max ,e = e e = 10 TeV, E min>p = j$m p v$ s = 14 keV, and E maXtP = e p = 
10 TeV. Similarly, 

n Hc fE m ^ e dEg(E,m e ,r e ,e e ) 
A^e = A e g (7 J 

Ue JeZZ'h' dEg(E,mne,Tue,en c ) 
= 1.0 x 10~ 7 cm" 3 GeV r ™\ (8) 

where E m - m: u e = ^m,He^f s = 57 keV and E maX: n e = eu e = 20 TeV. The parameters of the cosmic- 
ray proton and helium spectra are listed in Table 5, and the spectra are shown in Figure 9. These 
parameters and equation (2) yield a cosmic-ray injection efficiency r\ = 5 x 10 -4 if the magnetic 
field B = 40 fiG. This efficiency is assumed to be the same for all particle species. Our assumptions 
about the shapes of the cosmic-ray spectra and about the relative numbers of nonthermal protons 
and electrons yield a number density of protons that is 160 times larger than the number density 
of electrons at 1 GeV (i.e., dn c p jdE = 160 dnf/dE at 1 GeV; Fig. 9). Within the rather large 
uncertainties of our estimates, this ratio is consistent with the ratio observed at Earth (Meyer 
1969). 

The energies in cosmic-ray electrons, protons, and helium nuclei can be obtained from equa- 
tion (2) and the particle-spectra parameters listed in Table 5 and above: 

/*£/ m ax,e 

E e = V dEEA e g{E,m e ,T e ,e e ) 

" J - / min,e 

= 9 x 10 47 ergs, (9) 

/*_C/ m ax,p 

E p = V dEEA p g(E,m p ,T p ,e p ) 

" -C/ m in,p 

= 1 x 10 50 ergs, and (10) 

E-Rc = V dEEA He g(E,m.}ie,^}ie,eKe) 

= 8 x 10 48 ergs. (11) 
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The sum of these three energies is E CT = 1 x 10 50 ergs. As expected, the cosmic-ray energy is 
dominated by the energy in cosmic-ray protons. The estimates of the cosmic-ray and magnetic 
field energies are listed in Table 4 for the case in which B = 40 /xG. Dyer et al. (2001a) estimate 
the energy in only cosmic-ray electrons and obtain an energy that is about an order of magnitude 
larger than 9 x 10 47 ergs because they used a magnetic field of about 10 /iG instead of the field 
strength of 40 zzG used here. 

Drury et al. (1989) and Markiewicz et al. (1990) modeled the evolution of a supernova remnant 
using conditions that are quite similar to the conditions appropriate for SN 1006. These authors 
assume that 1 M Q of material is ejected with an initial kinetic energy of E$ = 1-2 x 10 51 ergs into 
an ambient medium that has a density uq = 0.1 cm~ 3 . The predicted results for the case in which 
the injection efficiency e = 10~ 3 (Fig. 3 of Drury et al. 1989 and Fig. 4 of Markiewicz et al. 1990) 
are consistent with our estimates of the distribution of energy between thermal particles, cosmic 
rays, and kinetic energy (Table 4). Note that the cosmic-ray injection efficiency parameter e of 
Drury et al. (1989) and Markiewicz et al. (1990), which should not be confused with our use of e 
for the exponential cutoff energy of the cosmic-ray spectra, is computed somewhat differently from 
the injection efficiency parameter r/ used here. If B = 40 /iG, ij = 5x 10~ 4 and the corresponding 
value of e = 1 x 10~ 3 . Therefore, our results and those of Drury et al. (1989) and Markiewicz et al. 
(1990) are consistent if the magnetic field is about 40 /iG. If the magnetic field strength is as small 
as 10 ifG, 7] = 2 x 10~ 3 , e = 4 x 10~ 3 , and the total energy in cosmic rays E cr = 4 x 10 50 ergs 
(i.e., the energy is expected to be more or less equally distributed between kinetic energy, thermal 
energy, and cosmic rays). For comparison, Figures 2 and 3 of Drury et al. (1989) suggest that 
almost all of the energy would be contained in the cosmic rays if e = 4 x 10~ 3 . Therefore, the 
results of Drury et al. (1989) are more compatible with a magnetic field of 40 /iG than a field of 
10 /iG. 

The estimates of the parameters of the cosmic-ray spectra of SN 1006 have important im- 
plications for Galactic cosmic-ray acceleration. It is generally accepted that Galactic cosmic rays 
are accelerated predominantly in the shocks of supernova remnants. If this hypothesis is true, an 
average supernova remnant should transfer a sufficient amount of energy to cosmic rays, produce 
a cosmic-ray proton spectrum that has an appropriate spectral index, and accelerate protons to 
energies high enough to explain the properties of the cosmic rays observed at Earth. Since the 
average flux of cosmic rays in the solar system has been constant within a factor of 2 over the last 
10 9 yr (Reedy, Arnold, & Lai 1983), the mean rate at which cosmic rays are energized should be 
approximately equal to the rate at which energy is lost as cosmic rays diffuse from the Galaxy. 
This rate (E loss = 0.3-1 x 10 41 ergs s _1 ; Parker 1969; Blandford & Eichler 1987; Drury et al. 1989; 
Lingenfelter 1992) and the assumption that one supernova remnant occurs every 30 yr suggest that 
an average supernova remnant produces 0.3-1 x 10 50 ergs of cosmic rays over the lifetime of the 
remnant (i.e., 3%-10% of the initial kinetic energy of the ejecta is transferred to cosmic rays). The 
estimate of the cosmic-ray energy in SN 1006 at the present is already consistent with this range if 
B < 100 /iG. Therefore, the total energy of the cosmic-ray protons in SN 1006 may be consistent 
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with the energy inferred for the sources of Galactic cosmic rays. 

The average spectral index of the relativistic cosmic-ray proton spectra produced by the ac- 
celerators of Galactic cosmic rays (T„ = 2.2) is inferred from the spectral index of the cosmic-ray 
proton spectrum observed at Earth (T p = 2.80 ±0.04; Asakimori et al. 1998) and the spectral steep- 
ening due to the energy-dependent escape of cosmic rays from the Galaxy (Ar = 0.6; Swordy et al. 
1990). The radio synchrotron flux of SN 1006 is produced by electrons that have energies of about 
1 GeV. At these energies the spectral index of the cosmic-ray electrons T e = 2.14 ± 0.12. Since 
relativistic protons and electrons are thought to be accelerated in a similar manner, the relativistic 
proton spectrum of SN 1006 may have a spectral index that is consistent with the index inferred 
for the accelerators of Galactic cosmic rays. 

The spectral index of the all-particle cosmic-ray spectrum observed at Earth is more or less 
the same up to an energy of about 3000 TeV. Above this energy, the spectrum gradually steepens 
to have an index of about 3. Since the observed cosmic rays that have energies below this "knee" 
feature are thought to be Galactic, the Galactic cosmic-ray accelerators are expected to be capable 
of accelerating particles to energies as high as 3000 TeV. If the mechanism responsible for the 
acceleration of Galactic cosmic rays depends on the magnetic rigidity of the particles (which is 
true for diffusive shock acceleration in supernova remnants), it may be the case that the cosmic- 
ray particles at 3000 TeV are principally iron and that protons are accelerated to energies of only 
about 100 TeV (Lagage & Cesarsky 1983). The estimated cutoff energy of the electron spectrum 
of SN 1006 (e e = 10 TeV), while uncertain, is lower than 100 TeV. Since relativistic electrons and 
protons that have the same energy have the same magnetic rigidity, the maximum energy of the 
protons is expected to be the same as the maximum energy of the electrons unless the maximum 
energy of the electrons is regulated by radiative losses. If the magnetic field B = 40 /iG, an electron 
with an energy E = e e = 10 TeV radiates half of its energy in about 400 yr. Since this time is less 
than the age of the remnant, the maximum energy of shock- accelerated electrons may be limited 
by synchrotron losses. In this case, the value of the energy e e represents a lower limit on the 
exponential cutoff energy of the proton spectrum because radiative losses are only important for 
electrons, not nuclei. Therefore, the maximum energy of the protons may be consistent with the 
expected value of e p = 100 TeV. However, if the magnetic field strength B = 10 /uG, the time 
required for a 20 TeV electron to radiate half of its energy (3000 yr) is significantly greater than 
the age of the remnant. If the magnetic field is this small, the maximum energy of the electrons 
in SN 1006 may be regulated by the free escape of the particles from the remnant (Reynolds 1996; 
Dyer et al. 2001a). In this case, the maximum energy of the protons is the same as the maximum 
energy of the electrons and is well below the expected energy of 100 TeV. 



5. Conclusion 

We present the results of a spectral analysis of RXTE PCA, ASCA SIS, and ROSAT PSPC 
data of the supernova remnant SN 1006. These data were fitted with several sets of thermal and 
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nonthermal X-ray emission models to characterize the global spectral properties of the remnant. 
The present work represents the first attempt to model both the thermal and nonthermal X-ray 
emission over the entire X-ray energy band from 0.12-17 keV. 

The best-fit model includes a nonequilibrium ionization component and a broken power law 
component, which are absorbed through an interstellar column density nu = 5.6 x 10 20 atoms cm~ 2 . 
The thermal component is described by an electron temperature kT e = 0.6 keV and an ionization 
timescale not = 9 x 10 9 cm~ 3 s. The large relative elemental abundance of silicon (10-18 times 
larger than the relative elemental abundance of silicon in the solar system) indicates that most of 
the X-ray-emitting silicon is reverse-shocked ejecta. The inferred mass of the shocked material is 
several times larger than the Chandrasekhar mass, which implies that most of the shocked material 
is interstellar. 

The nonthermal X-ray emission is well described by a broken power law component with a low- 
energy photon index Ti = 2.1 and a high-energy photon index 1^ = 3.0. Since 1^ is significantly 
larger than T\, the broken power law represents an approximation to a nonthermal spectrum that is 
steepening with increasing energy. This result supports previous claims that the nonthermal X-ray 
emission from SN 1006 is produced by synchrotron radiation from very high energy electrons. Using 
both the radio and X-ray synchrotron results, the spectrum of relativistic electrons in the remnant 
is inferred to have the form dn e /dE = A e E~ Te exp(— E/e e ), with A e = 2.4 x 10~ 9 cm -3 GeV r_1 , 
T e = 2.14 ± 0.12, and e e = 10 TeV. The values of A e and e e are based on the assumption that the 
mean strength of the magnetic field in the synchrotron-emitting region B = 40 fiG. This value was 
arbitrarily chosen for the present work because it is between the value obtained if the magnetic field 
is merely a compressed ambient field (10 uG) and the value obtained using the minimum-energy 
condition (100 uG). Other published estimates of the field strengths are near 10 /xG, but this value 
seems to be incompatible with the large, abrupt change in the radio flux at the bright rims, the 
detection of a significant radial component of the field, and the models of Drury et al. (1989). 

Since SN 1006 is expected to accelerate cosmic-ray nuclei as well as electrons, we have also 
estimated the parameters of the proton and helium spectra. The results suggest that the number 
density of protons is 160 times larger than the number density of electrons at 1 GeV. This ratio is 
comparable to the ratio observed at Earth. Integrating over the nonthermal particle spectra yields 
an estimate of the total energy in cosmic rays E CT = 10 50 ergs. This energy and the differential 
spectral index of the electrons at 1 GeV (r e = 2.14 ± 0.12) are consistent with the hypothesis that 
Galactic cosmic rays are accelerated predominantly in the shocks of supernova remnants. However, 
the maximum energy of cosmic-ray electrons in SN 1006 is well below 100 TeV (e e = 10 TeV). If the 
maximum energy of cosmic-ray protons e p = e e , the remnant does not accelerate particles to the 
energy of the "knee" at about 3000 TeV. Yet, if B = 40 uG, the maximum energy of the electrons 
may be limited by synchrotron losses and the maximum energy of nuclei might be significantly 
larger than 10 TeV. Therefore, SN 1006 is clearly a significant source of Galactic cosmic rays (at 
least cosmic-ray electrons), but the remnant may or may not be capable of accelerating particles 
to energies as high as the energy of the knee. Similar conclusions are reported for other young 
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Galactic shell-type remnants (Allen et al. 1997, 1999; Reynolds & Keohane 1999). 

Unlike the X-ray and radio synchrotron emission, TeV gamma-ray emission is only observed 
from the northeastern rim of the remnant. One possible explanation for this discrepancy is that 
the two rims have different electron spectra. The results of the present analysis are consistent with 
the idea that the electron spectra of the two rims are the same, but the possibility that the spectra 
are different cannot be excluded. Therefore, the nature of this puzzle remains a mystery. 
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Fig. 1.— A ROSAT PSPC image of SN 1006. North is up and east is to the left. The three 
rectangles are the regions used to extract the PSPC spectra for the northeastern rim (NE), center 
of the remnant (C), and southwestern rim (SW). The three ovals are the regions used to extract 
the ASCA SIS spectra. 
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Fig. 2. — RXTE PCA data of SN 1006. The top panel includes data for pointings at the the 
northeastern (NE) and southwestern (SW) rims of the remnant. The entire remnant is in the field 
of view of the PCA at both pointing positions. The histograms through the data show the results of 
the best-fit model. The spectral data and model of the southwestern pointing have been multiplied 
by a factor of 0.3 for clarity. The bottom panels show the ratios of the data to the model for both 
pointings. The model describes the spectra quite well. 



27 



1.000 



i> 0.100 

CD 
CO 

if) 



g 0.010 



0.001 




1.5 




1.0 


"0 


0.5 


"D 




o 


1.5 


<^_ 


1.0 


CC 


0.5 


Ctf 




Q 


1.5 




1.0 




0.5 



rf ,+++ 



^Aa#A#V. 



+ +■" i +i 



NE 



ASCA SIS 




, VW V ( ^^^^^ 



sw 



^ t+/+ XyVW¥^ 



0.5 



1 2 

Energy (keV) 



28 



Fig. 3. — ASCA SIS data of SN 1006. The top panel includes data for spatially separate regions 
along the northeastern (NE) and southwestern (SW) rims and in the center of the remnant (C). 
The histograms through the data show the results of the best-fit model. The spectral data and 
models of the southwestern and central regions have been multiplied by factors of 0.6 and 0.4, 
respectively, for clarity. The bottom panels show the ratios of the data to the model. In general, 
the model describes the spectra quite well. 
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Fig. 4. — ROSAT PSPC data of SN 1006. The top panel includes data for spatially separate regions 
in the center of the remnant (C) and along the southwestern (SW) and northeastern (NE) rims. 
The histograms through the data show the results of the best-fit model. The spectral data and 
models of the southwestern and northeastern regions have been multiplied by factors of 0.25 and 
0.10, respectively, for clarity. The bottom panels show the ratios of the data to the model. The 
model describes the spectra quite well. 
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Fig. 5. — 1,2, and 3 a confidence level contours for the plane of parameter space defined by two of 
the parameters of the nonequilibrium ionization thermal component. Here T s is the temperature 
associated with the forward shock and UqEq is the product of the square of the ambient density of 
hydrogen and the initial kinetic energy of the ejecta. The best-fit values of these parameters are 
indicated by the asterisk. From top to bottom, the three dashed lines are lines along which the 
ionization timescale not = 200, 500, and 1000 cm" 3 yr, respectively. From bottom to top, the five 
dotted curves are curves along which the characteristic electron temperature kT e = 0.3, 0.4, 0.5, 
0.6, and 0.7 keV, respectively. 
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Fig. 6. — 1, 2, and 3 a confidence level contours for the plane of parameter space denned by 
the photon indices of the broken power law. The best-fit values of the two indices are indicated 
by the asterisk. The dotted lines are lines along which the difference between the two indices 
Ar = T 2 - Ti = 0.55, 0.75, 0.95, 1.15, and 1.35. The nonthermal X-ray spectrum of SN 1006 
steepens with increasing energy since the line defined by Ar = is excluded at much more than 
the 3 a confidence level. This result supports the claim that the nonthermal X-ray emission is 
synchrotron radiation. 
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Fig. 7. — Radio to gamma-ray photon energy-flux spectrum of SN 1006. The data, which are 
labeled vertically, include the radio results of Kundu (1970), Milne (1971), Milne & Dickel (1975), 
Stephenson, Clark, & Crawford (1977), and Roger et al. (1988), the IRAS infrared upper limits 
of Arendt (1989), the ROSAT PSPC and RXTE PCA results of this paper, the EGRET gamma- 
ray upper limit of Hartman et al. (1999, Fig. 3), the gamma-ray results of the CANGAROO 
collaboration (Tanimori et al. 1998), and the gamma-ray upper limit of the JANZOS collaboration 
(Allen et al. 1995). The four model spectra are estimates of the photon energy fluxes produced 
by synchrotron radiation (S), inverse Compton scattering on the cosmic microwave background 
radiation (IC), the decay of neutral pions (7T°), and bremsstrahlung emission of the nonthermal 
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electrons (NB). The nonthermal X-ray spectrum is consistent with a model of the synchrotron 
spectrum, but the data are not consistent with the models of nonthermal bremsstrahlung and 
inverse Compton scattering. 
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Fig. 8. — Comparison of the nonthermal X-ray spectra of the northeastern and southwestern rims 
of SN 1006 in the ROSAT PSPC energy band. Here Tne — Tgw is the difference between the photon 
indices of the two rims and -Fsw/^ne is the ratio of the nonthermal fluxes at 1 keV. The best-fit 
values of the two parameters are indicated by the asterisk. The low-energy nonthermal spectra 
of the two rims are consistent with each other because the point defined by Tne — Tsw = and 
i^Sw/^NE = 1 is inside the 1 a confidence level contour. From the right to the left, the five dashed 
curves are curves along which Fsw/^ne = 0.5, 0.75, 1, 1.5, and 2.5, respectively, at 0.1 keV. 
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Fig. 9. — Estimates of the cosmic-ray electron, proton, and helium spectra of SN 1006. The low- 
and high-energy ends of the electron spectrum produce the observed thermal and nonthermal X-ray 
emission, respectively. The GeV electrons produce the observed radio emission. The ratio of the 
number densities of protons and electrons at 1 GeV is about 160, which is consistent with the 
ratio observed at Earth. The total cosmic-ray energy is dominated by the energy of the cosmic-ray 
protons. 
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Table 1. Results Obtained Using Two-Component Models 



n H 



Model 3 



(10 ° atoms cm 2 ) 



Component 1 
(Dimensionless or keV) 



Component 2 
(keV or cm~ 3 s 



) X 2 /v 



BPL+NEI 



5.60; 



-0.58 
-0.51 



Ti = 2.08" 

E b 

r 2 = 3.02 



0.11 

0.13 

i oc+0.18 

J-^-o.ie 

0.16 
0.14 



kT e = 0.58lg;§4 1109/729 
n t = 8.6t 7 o 7 7 x 10 9 



BPL+RS 


r oo+O.Vl 
°- oo -0.60 


ri = 2.06±&i£ 
E b = 1.90±g;i| 

r 2 = 3.1418:81 


A;T e — 0.162_ ; 009 


1202/729 


PL+NEI 


6.48^ 


T — 9 K1+U-04 


fcT e = 0.58 b 
n t = 8.6 x 10 9b 


1864/731 


Brem+NEI 


c i q+0.65 
D - id -0.57 


ft'-'ei — z - ly -0.10 


kT e2 = 0.58 b 
n t = 8.6 x 10 9b 


2338/731 



a These X-ray spectral models include a power law (PL), a broken power law (BPL), a thermal 
bremsstrahlung model (Brem), a thermal model based on the work of Raymond & Smith (1977; 
RS), or a nonequilibrium ionization model (NEI; Hamilton et al. 1983). 

b The 1 a errors of the characteristic electron temperatures, kT e , and the ionization time- 
scales, not, of the NEI model were not determined for these two sets of models, but they are 
expected to be similar to the uncertainties listed for the BPL+NEI model. 
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Table 2. Parameters of the Best-Fit Model 3 



Parameter 


Value 


90% Interval 15 




Interstellar absorption component 


liH (10 20 atoms cm -2 ) 


5.60 


5.04-6.24 


Broken power law component 



E b (kcV) 

r 2 



Flux(0.1-2 kcV) (10- 10 erg cm" 2 s" 1 ) 



2.08 
1.85 
3.02 
1.42 



1.94-2.20 
1.66-2.05 
2.86-3.19 

1.31-1.54 



Noncquilibrium ionization component 



T fs (K) 




10 8 


UqEq (ergs cm -6 ) 




10 5() 


Hc d 




1 


C d 




1 


N d 




1 







1.6 


Nc 




0.3 


Mg 




2.0 


Si° 




16 


s e 




16 


Ca c 




16 


Fe d 




1 


Ni d 




1 


Flux(0.1-2 keV) (10~ 10 erg cm- 2 


S" 1 ) 


2.53 



10' 



>-10 8 



10 49_ 10 50 



0.8-2.8 
0-0.8 
0-5.9 
10-18 
10-18 
10-18 



2.25-2.91 



a This model is the "BPL+NEI" model of Table 1. 

b The 90% confidence level intervals are computed using only the statistical uncertainties. 

c dF/dE = K(E/1 kcV)- r i for E < E b and dF/dE = K(E b /l kcV) r 2- r i (E/l kcV)- r 2 for 
E> E b . 

d The relative abundances of the elements helium, carbon, nitrogen, iron, and nickel are fixed to be 
the same as the relative solar abundances (Anders & Grevesse 1989) of these elements. 

e The relative abundances of the elements sulphur and calcium are fixed to be the same as the 
relative abundance of silicon. 
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Table 3. Parameters Inferred from the Best-Fit Thermal Model 



Parameter 


Value 


90% Interval* 


kT e (keV) 


0.58 


0.31-0.60 


Vf s (km s _1 ) 


2700 


1000-2900 


T{ S (pc) 


5.9 b 


5.6-16.3 


d (= r is /8) (kpc) 


1.4 b ' c 


1.3-3.7 


not (10 9 cm -3 s) 


8.6 


7.6-20.6 


t(yr) 


870 b 


760-6300 


n (= (nlEo/Eo) 1 / 2 ) (cm" 3 ) 


0.32 b 


0.10-0.37 


E (10 51 ergs) 


1.00 b 




M s (M ) 


9.7 b 


8.3-67 



a.r 



1 The 90% confidence level intervals are computed us- 
ing only the statistical uncertainties. 

b These parameters depend on Eq, which is assumed 
to be 10 51 ergs. The parameters rf s , d, and t 
scale as (-Eo/10 51 ergs) 1 ' 2 . The parameter hq scales 
as (E'o/lO 51 ergs) -1 ' 2 . The parameter M s scales as 
(-Eo/10 51 ergs). 

c The angular radius 8 = 15'. 
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Table 4. Inferred Distribution of Energy 





Value 


Parameter 


(ergs) 


E kin = \Mv 2 « ^M s vl 


4 x 10 50a '' 


EkT p = ^NpkTp « |7rm p nor|i;| 


4 x 10 50b 


£fcT e = f^e^Te « 27r^n r|fcT e 


2 x 10 49b 


-Eb 


3 x 10 47c 


-E'er 


1 x 10 50d 


Distribution of cosmic-ray energy E CT c ' d 


E e 


9 x 10 47 


E p 


1 x 10 50 


Ekc 


8 x 10 48 



a The kinetic energy of the unshocked ejecta 
is neglected. 

b These parameters scale as (Eq/10 51 ergs). 

c The magnetic field strength in SN 1006 is 
assumed to be 40 //G. 

d Within an accuracy of about 25%, these 
parameters scale as (13/40 /uG) -1 . 
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Table 5. Parameters of the Cosmic-Ray Spectra 





A 




e 


Particle 


(cm" 3 GeV 1 ^ 1 ) 


r 


(TeV) 


Electrons 


2.4 x 1(T 9 


2.14 ±0.12 


10 


Protons 


1.1 x 1(T 6 


2.14 


10 


Helium 


1.0 x 1(T 7 


2.14 


20 



Note. — The cosmic-ray injection efficiency rj = 5 x 
10~ 4 and is assumed to be independent of the particle 
species (see sec. 4.3). The ratio of the number density 
of cosmic-ray protons to the number density of cosmic- 
ray electrons is 160 at 1 GeV (Fig. 9). 



